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Abstract 

A soliton cellular automaton associated with crystals of symmetric 
tensor representations of the quantum affine algebra Ug(A^j) is intro- 
duced. It is a crystal theoretic formulation of the generalized box-ball 
system in which capacities of boxes and carriers are arbitrary and in- 
homogeneous. Scattering matrices of two solitons coincide with the 
combinatorial R matrices of U'JA^_^). A piecewise linear evolution 
equation of the automaton is identified with an ultradiscrete limit of 
the nonautonomous discrete KP equation. A class of N soliton solu- 
tions is obtained through the ultradiscretization of soliton solutions of 
the latter. 

1 Introduction 

The box-ball system invented by Takahashi and Satsuma |T5|] is an important 
example of soliton cellular automata. It is a discrete dynamical system in 
which finitely many balls move along the one dimensional array of boxes 



under a certain rule. Its integrability has been proved in ||TTMS|| by making 
a connection to the difference analogue of the Lotka-Volterra equation [H'T| 
through the limiting procedure called ultradiscrertization. 
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By now the original box-ball system has been generalized into several 
directions. First, one can introduce the balls distinguished by the index 
from the set {1,2,... , M}. Second, one lets the box at site n accommodate 
up to 9 n balls, where the capacity 9 n may depend on n. Third, one can 
introduce a carrier with capacity K t to redefine the time evolution at time 
t. The carrier comes from the left and proceeds to the right, picking up the 
balls in a box and dropping them into another under a certain rule. While 
it goes through the array of boxes, the successive loading-unloading process 
induces the motion of balls over the boxes hence the time evolution of the 
system. These generalizations of the Takahashi-Satsuma box-ball system are 
characterized by the parameters (M, 9 n , n t ). (n,t G Z play the role of space 
and time coordinates as in the figure in Section 2~3.) The original one [TS] 
corresponds to the choice (M,\/9 n ,Wn t ) = (1,1, oo). The case (M,\/9 n = 
1, V«t = oo) was introduced in |Tj and studied in [ [I'M S| | . Similarly, the cases 
(M = l,V0 n = 0, VVc = k) with k > 9 and (M,9 n ,Wn t = oo) were treated in 
|TM|1 and [ TTM]| , respectively. These works have been done mainly from the 



viewpoint of the ultradiscretization. 

The purpose of this paper is to study the general (M, 9 n , K t ) case. In Sec- 
tion |2] we formulate the corresponding generalization of the box-ball system 
in terms of the crystal theory |K|, |KMN1| , |KMN2|| . The latter is a represen- 
tation theory of quantum groups at q = 0. The unexpected link between 
the crystals and the box-ball systems has also been exploited in ||HKT , FQY 



through a crystal theoretic interpretation of the L-operator approach | |H1K 
The idea is to regard the box-ball system as a solvable vertex model |B|] at 
q = under a 'ferromagnetic' boundary condition. More concretely, the box- 
ball system corresponding to the data (M, 9 n , n t ) is naturally related to the 
U' q (A^I) vertex model at q = whose inhomogeneity in the quantum and 
auxiliary spaces is parametrized by 9 n 's and K t 's, respectively. 

Let B\ be the classical crystal of U' q (A^)) in the sense of [[KJV1JN1 



corre- 
sponding to the /-fold symmetric tensor representation of U q (AM)- Then the 
array of boxes and the ball configurations are identified with the elements 
from • • • ® Bg n ® Bg n+1 eg) • • • . The time evolution by the carrier with capac- 
ity Kt is realized as the action of the q = row transfer matrix acting on 
• • -®Bo n ®B$ n+1 ®- ■ • with the auxiliary space corresponding to B Kt . We call 

the resulting dynamical system the automaton. It is the most general 



one in the Aj$ case as far as the crystals for symmetric tensors are concerned. 
For generalizations to other root systems, see [IFTJ for a supersymmetric one 
and |[HKT|| for the non exceptional series other than A 

In Section |3| we introduce solitons and study the 2 soliton scattering. 
As in [FTKT, |FUY| we label the solitons in terms of the elements of the 



(i) 

M ■ 
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^g(^M-i)" ci y s t a l Bi, where I plays the role of the amplitude of a soliton. In 
the collisions of two solitons associated with B\ and B k , the scattering matrix 
is shown to coincide with the combinatorial R matrix giving the isomorphism 
Bi® Bk — Bk ® B\ of the [^(A^ f ' ) _ 1 )-crystals. These features are essentially 
the same with the \/9 n = 1 case ||TJN S| , [HK1| , |FU Y| ] . A new aspect here 
is that depending on the amplitudes /, k and the parameters 6 n , K t , smaller 
soliton can overtake the larger one. This is most transparently understood 
by viewing the scattering from the cross channel. By interchanging 9 n and 
Kt, one can swap the role of the space and time and thereby the boxes and 
carriers. Then the curious scattering mentioned above reduces to the 'usual' 
one in the cross channel where the larger soliton overtakes the smaller one. 



In Section |3.4| we also give a brief sketch of the conserved quantities of our 



automaton following [[FOY 



In Section [| we set up piecewise linear equation for the relevant combi- 



natorial R matrix ||NY|| and the resulting evolution equation for the A$ au- 



tomaton. Extending the earlier result ||1"1'1VI| ] , we identify the evolution equa- 
tion with an ultradiscrete limit of the nonautonomous discrete Kadomtsev- 
Petviashivili (ndKP) equation. A class of iV soliton solutions is obtained 
through an ultradiscretization of the r functions. As in the previous case 



TTM|| one needs to make a fine adjustment of the fermion momenta enter- 
ing the vacuum expectation value expression of the r functions. Each soliton 
in the automaton is obtained by letting M solitons in the ndKP merge to- 
gether in the ultradiscrete limit. 

Section [| is a summary. Appendices [A] and [E| contain the details of the 
proofs of Proposition and Theorem POL respectively. 



2 Automata from crystals 

2.1 ^(^)-crystals 

Let Bk be the classical crystal of U' q {A^) corresponding to the fc-fold symmet- 
ric tensor representation. As a set it consists of the single row semistandard 
tableaux of length k on letters {1,2,... , M + 1}: 



Bk 



{ 



mi ■ ■ ■ rrtk 



mi e {1, . . . , M + 1}, m a < • • ■ < m k }, 



where we have omitted the k — 1 vertical lines separating the entries. We 
also represent the elements by the multiplicities of their contents. Namely, 
b = I mi • • • mi, | G B k is also denoted by b = (xi,x 2 , ■ 



4{i 



mi 


■■m k 


mi 


= <}■ 



.xm+i) with Xi 
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Denote the Kashiwara operators of B k by fa and e, for i — 0, 1, . . . , M. 
The actions of e^, fa on B k are defined as follows: for b = (x±, x 2 , ■ • ■ , xm+i) £ 



e b = (xi - l,x 2 , ■ ■ ■ ,x M +i + 1), 

fob = (xi + l,z 2 , • • • ,x M +i - 1), 

ejfe = (a?!, • • • + - 1, • • • ,x M +i) fori = l, 

fib = (xi, • • • ,Xi - l,x i+1 + 1, • • • ,x M +i) fori = 1, 



,M, 
M. 



In the above, the right hand sides are to be understood as if they are not in 
B k . A crystal can be regarded as a colored oriented graph called a "crystal 
graph" by defining 



b' 



fib = b'. 



Thus for example B\ = {[T], . . . , M+l } has the crystal graph: 









M-l 


M 


M , 


M+l 







Setting Si{b) = max,{e{6 + | I > 0} and ^(6) = maxj{^6 + | I > 0} 
for b E B k , one has 

£ (&)=£i, £i(&) = x i+ i fori = 1, • • • , M, 
<Po(b)=x M +i, <Pi(b)=Xi fori = 1, • • • ,M. 

This data is necessary when we treat tensor products of the crystals. For 
two crystals B and B', the tensor product B <g> B' is defined. As a set, 

B <g> 5' = {6i <g> 6 2 I 6i e 5, 6 2 G 5'}. 

The actions of and fa are defined by 

ei(h <S> b 2 ) 

fi(bi®b 2 ) 

Here ® 6 and 6 ® are understood to be 0. For two crystals B and B', the 
tensor products B' ® B and B ® B' constructed as above are again crystals 
which are canonically isomorphic. The isomorphism R : B' ® B ^ B ® B' 



e-ibi <S> b 2 if > £j(6 2 ) 

&i (g) ei& 2 if < £i(& 2 ), 

/i6i ®_b 2 if y»i(6i) > £j(fc 2 ) 

h®fab 2 if y>i(6i) < £i(6 2 ). 



(2) 
(3) 
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is called the combinatorial R matrix |[KMN1| , [N Y|| . By the definition R 
commutes with f\, gj for any i — 0, 1, . . . , M. (More precisely one introduces 
affine crystals and the associated energy function, but in this paper we shall 
exclusively treat classical crystals and concern the energy function only in 
connection with the conserved quantities in Section |T 



Example 2.1. Figure [TJ and Figure [2] are the crystal graphs of U'JA^] 
crystals B 2 <8> B 1 and Bi ® B 2 , respectively. 




Figure 1: crystal graph of U'JA^ ^-crystal B 2 ® B\ 




Figure 2: crystal graph of [/'(A^-crystal i?i ® i?2 
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Example 2.2. Let B' = B 2 , B = B x of [/'(A^-crystals 





13 


® 


2 


i->|l|<g> 23 


ii) i? : 


23 


<8> 


2 


i->|3|®|22| 



These are obtained by comparing the crystal graphs in Example |2.1| . 

We write the highest weight element in B k with respect to U q (Am) as u k . 

k 

u k = \n...i\= (fc,o,--- ,o) eB k . (4) 



2.2 Isomorphism 

Here we give an explicit procedure to obtain the isomorphism R : B k ® B\ — ► 
5; ® -Bfc without drawing the whole crystal graphs of B k ® Bi and Bi® B k . 

Let b\ ® b 2 be an element in B k ® B t such as &i = (xi, . . . , xm+i) and 
^2 = {vii ■ ■ ■ i Um+i)- We represent b\ <8> 6 2 by the two column diagram. Each 
column has M + 1 rows, enumerated as 1 to M + 1 from the top to the 
bottom. We put %i (resp. yi) dots • in the z-th row of the left (resp. right) 
column. 



bi ® b 2 



• ■ ■ 


• 






• ■ ■ 


• 




• ■ ■ 


• 




1 



• ■ ■ 


• 


VI 




• ■ ■ 


• 


V2 






• ■ ■ 


• 
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Proposition 2.3. The rule to obtain the isomorphism R is as follows. 

(1) Assume k > I (resp. k < I). Pick any dot, say » a , in the right 
(resp. left) column and find its partner m' a in the left (resp. right) col- 
umn. The 9' a is chosen from the dots which are in the lowest (resp. high- 
est) row among all dots whose positions are higher (resp. lower) than 
that of » a . If there is no such dot, we return to the bottom (resp. top) 
and the partner m' a is chosen from the dots in the lowest (resp. highest) 
row among all dots. Connect » a and »' a by a line. We call the lines in 
the latter case winding and in the former case unwinding. 

(2) Repeat the procedure (1) for the remaining unconnected dots (I — 1)- 
times (resp. (k — 1) -times). 
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(3) The isomorphism R is obtained by sliding the remaining (k—l) (resp. (I— 
k) ) unpaired dots in the left (resp. right) column to the right (resp. left). 

The R obtained by this rule has the correct property as the isomorphism. 
This fact has been proved in section 3 of ||NY|| . We will write the relation 

R : u (8 v i— > v' <8 u' also as u <8> v ~ v' (8 u' . Obviously one has 



U k ®Ui~Ui® u k 



(5) 



for the element 



Example 2.4. Let M = 2, k 
by the following diagrams: 



2, I = 1. Example |2.2j (i),(ii) are obtained 




The line in (i) is unwinding and that in (ii) is winding. 

Suppose b®b' G B®B' is mapped to b'®b G B'®B under the isomorphism 
B <g> B' ~ B' ® B of L^(A^ } )-crystals. A Z- valued function H on B <g> B' is 
called an energy function if for any % and b®b' G B®B' such that ej (6(8)6') ^ 0, 
it satisfies 



Hfaib ® 6')) = F(6 ® 6') + 1 if i = 0, ^ (6) > £ (6'), ^o(6') > e (6), 
= H{b ® 6') - 1 if i = 0, <^o(6) < £o ( & '), ^(g/) < eo (6), 

= i/(6 (8 6') otherwise. 



(6) 



When we want to emphasize B®B', we write Hbb 1 for H. This definition of 
the energy function is due to (3. 4. e) of [N Y|, that is a generalization of the 
definition for B = B' case in |[KMN1|| . The energy function is unique up to 
additive constant, since B (8 B' is connected. By definition, HgB'ib ® 6') = 
B-B'B^b 1 <8> 6). Throughout this paper we normalize it as 



H Bl B k {ui®u k ) = 0, 



(7) 



irrespective of / < k or Z > k. Then it is the result of [|NY|1 that the energy 
function is (—1) times the number of unwinding lines in the sense of Example 



2.4 



With a successive application of i?'s, one interchanges the order of tensor 
product pairwise and obtains the isomorphism of B^ x (8 • • • <8 Bk n and Bk Pi <8 
• • • (8 Bk p for any permutation P. The compatibility of this construction 



is guaranteed by the Yang-Baxter equation obeyed by R. 
assertion follows easily from Proposition |2TB". 



The following 
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Proposition 2.5. Let k\, k<z, ■ ■ ■ G Z>i fte any sequence. Suppose b (g (g 
■ • — ci<E>- • -®Cn®b' is valid for some b' and q 's under the isomorphism 

Bi (g (gi • • • (g — -Bfci <g • • • <g -Bfc n ® -Sj- -For fl/ij/ b £ Bi, there exists n 
such that b' = U\ for all n> n . 

This property will be needed in constructing the automaton in Section 



2.3 Automaton 

Let • • • , 9-i, 6*o, 9\, ■ ■ ■ and ■ • • , Kq, K\, • • • be two sequences of positive 
integers. Denote the former indices by n, and the latter indices by t. Consider 
the 2D-lattice with n- and t- directions, 

n-direction • • • <g B 6n _ x <g> B 9n <g> -Be n+1 <£>••• 
i-direction •■■(g) 5 Kt _ 1 (g £? Kt (g B Kt+1 (g • • - . 

In terms of the box-ball systems, # n is the capacity of the n-th boxCand K t 
is the capacity of the t-th carrier. 

Draw t- constant lines horizontally, and n- constant lines vertically. Num- 
ber the former downward, and the latter to the right. At any horizontal or 
vertical line segment of the lattice, we inscribe an element of the crystals 
in the following way. At the point labeled by (t, n), we put b l n G B 6n on 
the upper line segment and v n G B Kt on the left line segment. Thus we 
have b n +l G Bg n on the lower line segment and v n+1 G B Kt on the right line 
segment. 



n— 1 n n+1 



t-1 



t+1 













b l 

V* n 


<4+l 






b t+l 















We impose the condition that they are related by the combinatorial R matrix, 

R : vi <g bl 



bit 1 ® v l 



n+1- 



In the following sections, we consider the time evolution of the system down- 
ward. In view of Proposition ETB we can and will exclusively consider the 
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case where for any t, b l n 7^ ug n only for finitely many n's and similarly for any 
n, v l n 7^ u Kt only for finitely many Vs. Sometimes we ignore i^'s and display 
the time evolution of the system only with the arrays 

... b°_ 2 b\ b° b° b° 2 ... 
... b\ 2 6i x b\ b\ b\ ... 
... bl 2 b\ b\ b\ b\ .... 

In short, the evolution of the array {b^} to is determined by 

B Kt ® (• • ■ ® B 9n ® B 9n+1 ® • • • ) ~ (• • • ® S 0n ® S 9n+1 ® • • • ) ® B Kt 
u Kt <g> (• • • <g> fo^ (8 g> • • • ) ~ (• • • <g> b l + l ® <g> • • • ) <g> 

under the successive applications of the combinatorial R matrices R : B Kt <g> 
A B 9j (g)B Kt . 

Setting p = • • • ® 6^ ® ® • • • , we denote the time evolution induced 
by u Kt as above by T Kt (p) = • • • ® ® 6^ ® • • • . Obviously the time 
evolutions are invertible, and due to (|^) they are commutative, 

T K T K , = T K ,T K . (9) 

In the rest of the paper, the 2 dimensional lattice on which the automaton 
is defined should be appropriately understood either as large but finite or 
formally infinite depending on the situation. 

The following observation will turn out useful in the sequel. 

Remark 2.6. Interchanging the role of 'space' and 'time', one can view 
as the evolution of the array ■ ■ ■ <8> t^+i <8> v * +1 <8> t^+i <S> • • • to the left as 

Te n (- ■ ■ ® <S ® <+i ® ^n+\ ® •••) = •■• ® 4 +1 ® < ® < X ® • • • ■ 

Example 2.7. Let M = 3, V0„ = 1 and V«t = 00. 

■ ■ ■ 111142113111111111111- ■ ■ 
■ ■ ■ 111111421311111111111- ■ ■ 
■ ■ ■ 111111114231111111111- ■ ■ 
■ ■ ■ 111111111124311111111- ■ ■ 

• - • 111111111112143111111- • - 

• • • 111111111111211431111- • • 

where i denotes This is a typical 2 soliton scattering. One can see that 
a soliton with amplitude I moves to the right with velocity / if separated 
sufficiently Hence the larger solitons overtake the smaller ones. (See Section 
|3~T| for the precise definition of the solitons and their amplitude.) 
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Example 2.8. Let M = 3, V0„ = 2 and Vac* = 1. 



14 11 11 33 11 11 11 11 11 11 11 
1 — | — 4 — | — 1 — | — 1 — |— 3 — | — 1 — | — 1 — | — 1 — | — 1 — | — 1 — | — 1 — | — 1 

11 14 11 13 13 11 11 11 11 11 11 
1 — | — 1 — | — 4 — | — 1 — [—3 — [— 1 — | — 1 — | — 1 — | — 1 — | — 1 — | — 1 — | — 1 

11 11 14 11 33 11 11 11 11 11 11 
1 — ^ — 1 — | — 1 — | — 4 — | — 1 — |— 3 — | — 1 — ^ — 1 — ^ — 1 — | — 1 — | — 1 — | — 1 

11 11 11 14 13 13 11 11 11 11 11 
1 — | — 1 — ^ — 1 — | — 1 — | — 4 — [— 1 — [—3 — | — 1 — | — 1 — ^ — 1 — | — 1 — | — 1 

11 11 11 11 34 11 13 11 11 11 11 
1 — |— 1 —\— 1 —\— 1 —\— 1 — |— 3 —\— 1 — |— 3 —\— 1 —\— 1 —\— 1 —\— 1 

11 11 11 11 14 13 11 13 11 11 11 
1 — ^ — 1 — | — 1 — | — 1 — | — 1 — | — 4 — [— 1 — ^ — 1 — ^3 — | — 1 — | — 1 — | — 1 

11 11 11 11 11 34 11 11 13 11 11 
1 — ^ — 1 — | — 1 — | — 1 — | — 1 — | — 1 — ] — 3 — ^ — 1 — ^ — 1 — [—3 — | — 1 — | — 1 

11 11 11 11 11 14 13 11 11 13 11 
1 — | — 1 — ^ — 1 — | — 1 — | — 1 — | — 1 — | — 4 — [— 1 — | — 1 — ^ — 1 — [—3 — | — 1 

11 11 11 11 11 11 34 11 11 11 13 
Here i and ij denote [T| and ij , respectively. We have depicted the both 



variables {b l n } and {v^}. This time 14 on the top left is the smaller soliton 
and 33 or 13® 13 is the larger soliton. Thus in terms of the {b^} variable, the 
smaller one overtakes the larger one as we go down the figure ending with the 
solitons 34 and 13. This is an opposite feature from the previous example. 
However in the space-time interchanged picture (Remark |2^), it reduces to 
the situation similar to Example |2.7| . Namely, in terms of the {v^} variable, 
the larger soliton overtakes the small one as ... 43 ... 3 .. . — > ... 4 ... 33 ... , 
as we trace the diagram from the right to the left. 



2.4 Equivalence with box-ball systems 

Our automaton can be viewed as a generalized box-ball system. One 
interprets the letter 1 in the tableaux as an empty space and the other letters 
2<z<M + las the balls with index M + 2 — i. The element b l n signifies 
the balls contained in the n th box with capacity 8 n at time t. Similarly 
v l n stands for the carrier with capacity K t . Then ([8]) tells that through the 
loading-unloading process, the box and the carrier change into and v l n+l) 
respectively. Sending the carrier through to the left, one has the time evolu- 
tion of the box-ball state •■■®b t n ® ® ■ • • into • • • <g> <g> ® • • • . 
For a concrete rule describing (|8]) in terms of the box-ball terminology, see 
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BBS scattering rule in JTNSJ. The relation (|8]) will also be expressed as a 
piecewise linear equation in Proposition [4.1| . 

When V/% = oo we claim that the evolution of {b^} in our automaton 
is equivalent to the box-ball system studied in [[TTM]| under the above stated 
translation. In the latter the one-dimensional array of boxes with capacities 
... , n _i, 9 n , 9 n+ i, . . . accommodate the balls with an index from the set 
{1, . . . , My. The dynamics of the balls in each time step is governed by the 
rule ITTMfl : 



1. Move every ball only once. 

2. Move the leftmost ball with index 1 to the nearest right box with space. 

3. Move the leftmost ball with index 1 among the rest to its nearest right 
box with space. 

4. Repeat this procedure until all of the balls with index 1 are moved. 

5. Do the same procedure 2 — 4 for the balls with index 2. 

6. Repeat this procedure successively until all of the balls with index M 
are moved. 

If the ball with some index is absent, one just proceeds to those with the next 
index. A box with space means the one that contains strictly fewer balls than 
its capacity. If a box contains more than one balls with the same index and 
they are not yet moved at an instant during the procedure, one may pick 
any one of them when looking for the leftmost one. The equivalence to our 
automaton with V/% —>■ oo is shown by the fact that the both lead to the 



same evolution equation, which is given from Proposition |4.1| under the said 
limit. 

The above rule tells that the time evolution in our automaton admits 
the factorization: 

T oc = f M ---f 2 f 1 , (10) 

where Tj moves the balls with index j only, and we identify the left hand 
side with the corresponding operator acting on the box-ball systems. 

For a later convenience we introduce the canonical system following |[TNS|| . 
We keep assuming \/n t = oo and stay in the description in terms of the box- 
ball system rather than crystals until the end of this subsection. Thus we 
identify b G Bg with the capacity 6 box containing the balls as specified 
before. Suppose a state p = ■ ■ ■ <S> b n <g> b n+ \ (g) • • ■ contains J balls in total. 
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Then the action of Tm ■ ■ ■ T 2 T\ consists of J steps, each of which is to move 
a certain ball. To a ball to be moved in the j th step (1 < j < J), we 
assign a signature j. The assignment is unique up to the trivial freedom 
among the commonly indexed balls within the same boxes. Let c(p) be the 
ball configurations obtained from p just by regarding the signatures as new 
indices. It consists of the same array of the boxes and J balls as before but 
with the new distinct index from 1 to J. One can still let c(p) evolve under 
the previously stated rule 1 — 6 by replacing M therein with J. The resulting 
new box-ball system is called the canonical system. By a close inspection of 
the rule 1-6, it is not difficult to confirm the commutativity: 

c(f M ■ ■■f 2 f 1 {p)) = fj ■ ■■f 2 f 1 (c(p)). (11) 

In this sense the canonical system essentially grasps the time development 
pattern of the original one. This fact, firstly recognized in |[TNS|| for V# n = 1, 
will be utilized in Appendix |X[ 

3 Combinatorial R matrix as scattering ma- 
trix of ultra-discrete solitons 



Here we prove Theorem |3.10| , which identifies the scattering matrix of the 



ultra-discrete solitons with the combinatorial R matrix of U 'JA^p 



3.1 Solitons 

Let B' k be the classical crystal of U'(A$_^) corresponding to the fc-fold sym- 
metric tensor representation: 

B' k = {\m 1 ---m k \ | rrii G {1, ... , M}, m x < ■ ■ ■ < m k }. 



Denote the Kashiwara operators of B' k by f[ and e- for % = 0, 1, . . . ,M—1. For 
distinction, from now on we use the notation _B fc ,/j,ej for U'AA^j ')-crystals 
and B' k ,][,e! i for ^(A^J-crystals. Let R and R' be the combinatorial 
R matrices for U' q (A^) and U' q {A^_^)^ respectively. Thus Rfi = fcR and 
R'fl = flR' hold when they act on the tensor product of two crystals, and 
similarly for e^, e! v (We will specify the crystals that they act each time.) 

Remark 3.1. When M = 1 we still define B' k as above, which is the set with 



the single element u k = 1 ... 1 . We further understand that the 'V'(A 



combinatorial R matrix R' : B[ ® B' k — ■> B' k ® B[ is given by R'(ui <8> u k ) 



12 



For each k E Z>i define a map i k by 









(Si)®* 






m l ---m k 




m k + l 


® • • • <8> 


TOi + 1 



Let fci, . . . , fc/v G Z>i and L , . 
2^ we further introduce a map 



Ljv G Z>o for some N E Z>i. In terms of 



(£o,— j-E-at) . tdI 
l k lt ...,k N ■ -"fei 



si 



fejV 



^JV+fciH ^ fe iV 



by 



HP" ® Zfa (&i) ® 0^ ® Zfe (62) ® • • • ® Hf ^ 1 ® (M ® E 



In particular 
<jfc denote the map 



,(0,0) 



. The map zL 



(Lo,... ,L N ) 



ft: (Si)®* 



is injective. For each k E Z>i let 



Bi 



mi 



® •• • ® 







mi . 


•< 



where 1 < m' x < ■ • • < < M + 1 are just the re-ordering of mi, . . . , m k 
into the weakly increasing order. We assume that £ := J2 n n is sufficiently 
large. We set 



e 



Ba 



Ba 



For non-negative integers L , . . . ,Ln such that £ = L 

(L ,...,L ' 

t fel,... ,fcjv 



+ /cat, denote by 4f°'' fct^ ^he composition o z' L " /; x 



N 



(12) 
fei + 



fei,...,fcjv ' X ' e '' 



fc£° = ^ ® • • • ® Sl w 



^(L ,... ,L W ) 
fel )•-. ,fejy 



(13) 



Suppose that the image is obtained from the element ■ • • ® ug n ® ug n+1 <8> ■ ■ • 
by replacing only the isolated segments u$ n , ®Ug n . +1 ®- ■ - ®uq , (m < with 
some ba ® ■ • • <8> fee , G -Be,, <8> ■ ■ ■ <S> Bg , for 1 < % < N. Assume further that 

l % 

the interval is sufficiently large, namely, n, — n'^ ^> max(fei, . . . ./cjv) for any 
2 < % < N. In such a case we call the image of ([13|) an asymptotic N soliton 
state. Each soliton is essentially associated with an element in B' k , and we 
call k the amplitude of the corresponding soliton. States obtained from an 
asymptotic N soliton state under arbitrary time evolutions T K ■ ■ • T K i will be 
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called iV soliton states. This definition will naturally be justified from the 
consideration on the conserved quantities in Section p.4| . Note that 1^°' 
is not injective since is not. Consequently, the result of application of 
L ki' "k N N * s n °t necessarily an 'overall translation' of ( p~3|) in a naive sense 



even when L\ — Li is i-independent for i < N or i > 0. See Example [Ol 
below. 



First we consider N — 1 case. As it turns out in Proposition |3.3| , there 
is no distinction between an asymptotic 1 soliton state and a 1 soliton state. 
Moreover one can check that the definition of the 1 soliton state here agrees 
with the 1 soliton solution that will be given later in fl5T|). Given a 1 soliton 

state 



»n + l 



p = ■ ■ ■ <g> 6 n -i ®K® K+i ® • • • e • • • <8> -B 6)n _ 1 ® <8) B e 
one can unambiguously specify integers n,k(> l),s,t by the conditions: 
bj = Uq . if j < n or j > n + k, 



1 • • • lmi • • • m t 1 < t < n , 2 < mi < • • • < m t < M + 1, 
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1 • • • lmj • • • m' s < s < 6» n+fc - 1, 2 < mi < • • • < m' s < M + 1, 



b n+ i, . . . , do not contain 1 in their tableaux. 

Note that 'if in the first condition is not 'only if in that b n+ k = u e n+k 
is allowed as s = 0. The amplitude of the soliton according to the above 
definition equals t + n+ \ + • • • + n +fc-i + s - We set 

X (P) = i ~ y ^ = 1 + 0n+1 H h ^n+fc-l 

and call x(p) the coordinate of the soliton. should not be confused with 
the amplitude of the soliton. 

Example 3.2. Consider B 9l ® • • • ® B 06 with 0! = 3 = 1, 2 = 4 = 6 = 2 
and 5 = 3, hence £ = 11. 

(i) Take b = \\\eB[. Then ^"'^O) with Lq + L x = 10 are examples of 1 
soliton states with amplitude 1. One has ^'^(b) = if' 8 \b), ^'^(b) = ^'^(b) 
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and 4 M) (6) 



,(7,3) 



(b) 



,(8,2) 



(6). For L < 8 they look as 



P 



n n + k x(p) y(p) 



,(o,io) 
n 

,(2,8) 
,(3,7) 
,(5,5) 
,(8,2) 



iti o n^i <g) nn <g» pt] ( 

■ [T| (8) [lT| (8) |T| ® [12I ( 

[]]<8>[]]]®[]]<g>[T[]< 



111 



111 



111 



111 



112 



2 
3 
4 
5 

6 




2 
3 
5 



where we have also listed n,n + k, x(p) and yip). 

11223 I G B' 5 . Then i[ L °' Ll \b) with L + Lj 



(ii) Take c 



6 are examples of 



1 soliton states with amplitude 5. For Lq < 5 they look as 



p 














r? 


n + k 


x[p) 


y(p) 


4 0,6) ( C ) 


= 0^ 


S> |_33j (8) 


2 5 


D 12 




111 


g>0O 


1 


4 





4, 


4 1,5) (c) 


= 0^ 




3 £ 


9 22 




111 


»0D 


2 


5 


1 


5. 


4 2 ' 4) ( C ) 


= 0^ 




3]<E 


D 23 




112 




2 


5 


2 


4, 


4 3,3) ( C ) 


= 0^ 




3§ 


D | 33 | ® 


122 


^ 00 


3 


5 


3 


3, 


4 4 ' 2) ( C ) 


= 0^ 




05 




223 


» EH 


4 


6 


4 


5, 


4 M) ( C ) 


= 0^ 




05 


?[l4l® 


233 


g>pL2] 


4 


6 


5 


4. 



In Section |3.3| we will make use of 

Proposition 3.3. Let p = i[ L °' Ll \b) be the 1 soliton of amplitude I asso- 
ciated with b G B[. Then its time evolution T K {p) is again 1 soliton and 



expressible as T K (p) 



,Wi) 



(b) for some L' ,L[ (L' + L[ = L + L\) but 



with the same b G B[. The difference of their coordinates (velocity under T K ) 
is given by 



x{T K (p)) - x(p) 



K 



min(K, I) + max( 



'n+k 



k < y(p), 
1,0) n>yip). 



The proof is done by a cumbersome case study. When ^J9 n = 1, the above 
result simplifies to x(T K (jp)) — x(p) = min(fc, I) in agreement with [[FOY|] . In 
general, the velocity varies locally depending on the data {6 n }. In Example 

O (i) one has T K (4°' 10) (6)) = if' 8 \b),T K iif' 8 \b)) = if' 7 \b),T K (^ 7 \b)) = 



Li (b), T K (i^ ' o; (6)) = L ( i' z> ib) for any k > 1. Similarly in (ii) one has 



,(5,5), 



,(8,2) 



TJl 



(0,6), 



c)) 



(k',6-k') 



(c) for any k > 1, where k' = rnin(«, 5). These re- 



sults agree with Proposition EO 
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Let l 



Cl 



cjy) (cj G -£>£,.) be an asymptotic iV soliton state 



and 



• • • ® 61 ® 6* 



n+l 



-<■ Hi- 1 K t -1 



(Lq,— ,Ljv) i 



Cl 



Cat, 



be its time evolution. Assume that the solitons are enough separated without 
an interaction throughout the time interval in consideration. Let {v^} be the 
associated variables on the vertical edges as in (^). Then in the space-time 
interchanged picture, the state • • • Cg) v * +1 ® v * ® • • • is also an asymptotic iV 
soliton state associated with the same ci <8> • ■ ■ <S> cjv- Namely, 



KOI 



(L' ,-,L' N ) 



ki,... ,k N 



Cl 



for some L' r 



o- 

(• • • € s, 



, L N . 



Kt+l 



Here 



■ ■ ® C N ) 

<g> B Kt+1 ®B Kt ®- 



is an analogue of 6 in (0), and we have set M. = ^2t K t- The figure in 
Example |2.8| will be of help to understand this fact. In a sense one can employ 
either picture to describe the scattering process. Indeed our discussion in the 



end of Section 3.3 will rely on this observation. 



3.2 Scattering of 2 solitons; a typical case 



Our aim here is to show Theorem 3.S which is valid in the 'typical' situation 



(19). 



Lemma 3.4. For each i — 1, . . . , M — 1, we have a commutative diagram: 



Bl 



B' k □ {0} 



e-i+i 



[B,r k □ {o}, 



where = 0. The same relation holds also between f- and /, 



i+i- 



Combining Lemma |3.4| with the realization of Bq in Bf as a U q (AM) 
crystal (cf. |[KN|| ) , one can derive the following lemmas. 



Lemma 3.5. In the diagram 



B> ® B' k 



,( L 0' L 1> L 2) 



B dn ® B 9n+1 ® 



e i+ i 



,{ L 0> L 1> L 2) 

{B[ ® B' k ) U {0} J* > 



B 6n ®Be, l+1 ®---)U{0}, 
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suppose that the image of L lk ' 1: 2 is an asymptotic 2 soliton state. Then 
the diagram is commutative for any i — 1, . . . , M — 1. The same relation 
holds also between f[ and f i+ i . 

Actually, the commutativity of the above diagram holds under a milder 
condition than being an asymptotic 2 soliton state. 

Lemma 3.6. Letpx, . . . ,p m be the subsequence of ax, ... , (a n G Bg n ) con- 
sisting of all the elements such that a n ^ ug n . Assume the same relation 
between p[, . . . ,p' m and a[, . . . , a' L . Then for any t, t' G Z> and k G Z> 1? 
the two relations 

fi+i(Pi ® • • • ® Pm) = p[ ® • • • ®p' m , 

f i+ i{uf ® ai ® ■ ■ ■ ® a L ® wf) = iz®* ® ai ® • • • ® a' L ® 

are equivalent for each i — 1, . . . , M — 1. T7ie equivalence persists even when 
the right hand sides are both 0. Tae same is true also for ej+i. 



Proposition 3.7. Suppose an asymptotic two soliton state has evolved into 
another as 

n (it k oMM \b ® C )) = if'r L '^\d ® y) (M) 



for some K,t,Li,L[ > 0,6,6' G -B z ' and c, c' G Tnen ( fl^[ ) a/so va/id 
under the replacement of ' b ® c (Vesp. c' (g> 6'j 6y //(6 (g) c) (Vesp. f[{d ® b') ) 
for any 2 = 1,... , M — 1 st/cn iaai <g> c) 7^ 0. 

Proof. (0) is equivalent to 

«f ® $ MM \b ® c) ~ 4 i r i ' l,L ' ) (c' ® &') ® «f 



Apply to the both sides. Due to Lemmas |3.5| and |3.6|, the result becomes 



"... ' ■» C ^ ® c )) * 4? 1 2 (/i( c ' ® &/ )) ® «f • 



®* ^>> . ( L 0,il,i2) ,0. „\\ _, ,( L 0> L 'l< L 2) { f>{l 

□ 



Proposition 3.8. LeH > and assume that l^' u 2 (61(8)62) an asymp- 
totic 2 soliton state with 

6i = (Z,0,...,0)eB,', 6 2 = (h,k-h,0,... ,0) G 5£ (15) 
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with < h < k in the notation of (|I[). Assume further that I > 9 n for all but 
finitely many n 's. Then if k ^> I, there exists t > such that the result of 
the time evolution T* also becomes the asymptotic 2 soliton state as 

Tl (ti° M ' L2 \h ® 62)) = if L ' l ' L ' 2 \c 2 ® Cl ), (16) 

where Ci,c 2 are given by 

c 2 = (fc,0,... ,0) e B' k , ci = (l-k + h,k-h,0,... ,0) e B[. (17) 

The proof is given in Appendix |A|. In fact both 61 ® b 2 and C2 <S> ci 
are U q (AM-i) highest element, i.e., e-(&i ® 62) = e-(c2 <8> ci) = for all 
1 < % < M — 1. Combining this property with the conservation of weights 
(number of the letters) and the soliton content (cf. Section |3.4| ), one can 
argue that the outgoing state should necessarily correspond to c 2 ® ci if it 
is an asymptotic 2 soliton state at all. However, to establish the separation 
into 2 solitons asymptotically is not a trivial task for inhomogeneous ^„'s 
only bounded by the condition I > 6 n for all but finitely many n's. So far we 
have not managed it without recourse to the actual 2 soliton solution as in 
Appendix 0. 

As a C/ g (Ajvf_i)-crystal, the C/'(A^ ) _ 1 )-crystal B[®B' k decomposes into the 
connected components. Each component is parametrized with the U q {A.M-x) 
highest elements 61 ® b 2 (|T5|), and is generated by applying f[ operators 
(1 < i < M — 1) to it. The decomposition of the same pattern takes place 
also for B' k <S> B[ according to the highest elements c 2 <8> C\. Combining this 
fact with Propositions |3.7| and |3.8| , we conclude that there exists a map 5" 
(S matrix) uniquely defined by 

S' : B[ ®B' k ^ B' k ® B[ 

Tl ($ MM \* ® c )) = iT'^iS'ib <g> c)), (18) 
under the condition 



n ^> I > k, I > 9 n for all but finitely many n's. (19) 
It describes the 2 soliton scattering. 

,.. /in) ;. - ... r>' — c .... , ; \ -i r , 



Theorem 3.9. Under the assumption (\TQ), we have R' = S' on the U'JA^ 



crystal B[ ® B' k . 



Proof. By the definition and Proposition |3.7| , the both R' and S' commute 
with f[ for any 1 < i < M — 1. Moreover, for any U q (A M _i) highest elements 
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bi ( 
R'( 



n 



) 62 given by (|T5D, their actions are the same, i.e., S' 
~>i <S> b 2 ), where the latter can be verified from Proposition EO 



C 2 <X> Ci 



□ 



Thus in the situation ( [TT?D the larger soliton overtakes the smaller soli- 
ton and the scattering matrix coincides with the combinatorial R matrix of 
L^A^J -crystal. For instance Example |2j] tells that 

->pn®[23]. 



5': 13 



This agrees with Example |2.2| (i). 



3.3 Scattering of 2 solitons; general case 

First let us consider the homogeneous case V0 n = 0, Wet = «. Fix positive 
integers I > k. We study the scattering of 2 solitons in • ■ ■ ® Bg <g> Bg <g> ■ ■ • 
with amplitudes I and under the time evolution T*. The qualitative feature 
of the scattering depends on the cases: 

(i) I > k > max(#, k) 

(ii) min(^, k) > I > k 

(iii) l>K>k>6 

(iv) K>l>k>9 

(v) Z > k> 6>k 

(vi) K>l>9>k 

(vii) I > 9 > k > k, 

(viii) 9 > I > k > k, 

(ix) I > 9 > k > k, 

(x) 9 > I > k > k. 

Here the classification has been done so that 

{(i) n (ii)} II {(in) U (iv) U (v) U (vi)} H {(vii) U (viii) U (ix) U (x)} . 

For example (iii) and (iv) share I = k > k > 9 case. However the three groups 
are mutually disjoint and correspond to distinct features of the scattering as 
we will see below. The vi and Vk are the velocities of the solitons with 
amplitude I and k, respectively. For each soliton it has been calculated by 
using Proposition [3.3| by assuming no effect from the other soliton. In (vi) 
and (x) we have excluded I = 9 and I = k, respectively since they both lead 
to v 1 = v k = 9 hence no scattering. By the same reason the cases (i) and 
(ii) are out of question. Via the space-time interchange 9 <-> n, the cases 





Vl = v k 


= K 




Vl = v k 


= 9 


Vl 


= K > V k 


= k 


Vl 


= l> v k 


= k 


Vl 


= K > V k 


= 9 


Vl 


= I > v k 


= 9 
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(vii) ,(viii) , (ix) and (x) are mapped to (iii), (iv), (v) and (vi), respectively. 
(See the argument before Theorem |3.10| on the velocities in the cases (vii) 

- (x).) Thus we are left with (iii)-(vi), where / > 9 and v\ > Vk are always 
valid. Following [ FOY|1 , we utilize the commutativity (0) and consider the 2 
soliton scattering under T* as 

rjit rji — t' rjit rjit' 

ft OO K OO ' 

The scattering are thus divided into three stages. In the first stage, we 
let solitons evolve under for sufficiently large t' . Since / > 9 matches 
the condition fllTf), Theorem |3l^ tells that the larger soliton overtakes the 
smaller one with the scattering rule described by S' = R'. In the second 
stage corresponding to T*, the larger soliton goes further ahead than the 
smaller one with no interaction because of vi > Vk- Therefore in the last 
stage , the two remain isolated even though they are drawn back and 
get relatively closer. Thus we conclude that in all the cases (iii)-(vi), the 
qualitative feature is the same as the one in Theorem |3.9| . Namely, the 
larger soliton overtakes the smaller one and the scattering rule is given by 
the combinatorial R matrix R' : B[®B' k — > B' k ®B[. Through the space-time 
interchange argument, this implies the opposite feature of scattering in the 
cases (vii)-(x). Namely, the smaller one overtakes the larger one with the 
scattering rule given by the combinatorial R matrix R' : B' k ® B[ — > B[ ® B' k . 

We note that in the cases (vii) - (x) , one does not necessarily have vi < v k 
at any time. It actually depends on whether k ^ y(p) when Proposition [T3] 
is applied. Nevertheless v\ < Vk should be valid "on average" and the above 
feature of the scattering should hold due to the reduction to the cases (iii) 

- (vi) where the strict inequality v\ > Vk is always valid. To summarize, we 
have shown 

Theorem 3.10. Let I > k be the amplitude of 2 solitons in ■ ■ -^Bq^Bq®- ■ ■ . 
Under the time evolution T K , the scattering matrix of the collision (if any) in 
the sense of or ftltyi^k is given by S' = R' , where R' is the combinatorial 
R matrix of the U'^A^^)- crystals for 

(7) B[ ® B' k ^ B' k B[ if min(/, k) > max(k, 9), 
(II) B' k ® B[ A B\ ® B' k if min(/, 9) > max(k, k), 
(III) no scattering (same velocity) otherwise. 

Example [2.8| corresponds to the choice l = 9 = 2,k = n = l, hence to (II) 
in the theorem. The scattering matrix is read off the figure, 

S" : m O [221 [231 1^" . 
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This agrees with the inverse of the R matrix in Example |2.2| (ii). 

Let us comment on the inhomogeneous case where (9 n 's and K t 's actually 
depend on the indices. In view of flI5|), the qualitative feature of the scattering 
remains the same as Theorem |3.10| even if we slightly relax the conditions 
therein. For instance the larger soliton still overtakes the smaller one with 
the rule S' = R' if min(7, K t ) > max(&;, 9 n ) holds for almost all n and t that 
are relevant during the scattering in question. In such cases we expect that 
the asymptotic N soliton state in the sense of Section |0| undergoes the 
scattering which are essentially factorized into the two-body ones studied 
here. On the other hand, if 6* n 's and K t 's are not bounded by the condition 
as above and indeed far from being homogeneous, even 2 solitons can collide 
many times in general depending on the local velocities. In such a case we 
do not have a simple picture of the scattering. 

Example 3.11. Let M = 3. 
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where • denotes (g>, and 14 for example does 14 6 B 2 . Not only # ra 's but also Kt 
are inhomogeneous here so that the relevant time evolutions are T§ for the process 
0^1^2^3^4, whereas they are T 2 for 4 -> 5 -> 6 -> 7 -> 8 -> 9. This 
is an example of the double scattering of 2 solitons caused by the inhomogeneity. 
The larger soliton once overtakes the smaller one, but after the collision it gets 
slower due to the environmental change and is eventually passed by the smaller one 



again. This is easily understood from the classification (I)— (III) in Theorem 3.1C 
for the homogeneous case. In the first stage we have I = 4, k = 5, k = 2, 9 n < 3 so 
that the larger soliton overtakes the smaller as in (I). On the other hand we have 
k = 2, 9 n > 3 in the second stage hence the smaller one passes the larger one as in 
(II). Following the time evolution downward, one finds the scattering matrices for 
the successive collisions: 



1223 


«> 


13 


i — > 


23 


® 


1123 


i — > 


1223 


® 


13 



in terms of the soliton labels with the U'JA^ )-ciystal elements. They agree with 



the combinatorial R matrices B' 4 (£i B' 2 — B' 2 (?) B' 4 calculated from Proposition 2.2 
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3.4 Conserved quantities 



Let us give a class of conserved quantities in the A*) automaton. Since our 



construction here is based on ||FOY|| and the result is quite parallel, we will 



only present a brief sketch. Given an automaton state p 
■ ■ ■ (b n = u 9n for |n| > 1), let 



u K <g> p ~ • • • ® b' n _ x ®b' n ®v n 
for some b\ G Bq. and v n G B K . Set 

E K (p) = - H B K B Bn+1 (v„ ® b n+1 



+2 



which is well defined owing to the normalization (|7|). By the same argument 
IFDYI1 



as m 



we get 



E K (T K >(p)) = E K {p) for any k, k'. 



Thus E K (p), k G Z>i form a family of conserved quantities. If p is an 
asymptotic iV soliton state in the sense of Section j3J] , it is straightforward 
to derive 



E^p) =J^min(Z,K)iV I . 



(20) 



where iVj is the number of solitons with amplitude /. Therefore if a state with 
the soliton content {N{\ scatter into another state with the content {iV/}, 
Ni = N[ must be valid for any / due to the conservation of all E K 's. In both 
Example |2.7| and |2.8| , we have Ei = 2, Ei = 3 for I > 3, in agreement with 
N x = iV 2 = l,JVj = for Z > 3. In Example |T|, we have E 1 = 2,E 2 = 
4, E% = 5 and = 6 for Z > 4, in agreement with A^; = 812 + ^4. 

When \/8 n = 1, (|20f) is obtained in proposition 4.4 in [ FQY . An equivalent 



family of the conserved quantities has also been given in 

Another conserved quantity is the semistandard Young tableau, which can 
be constructed as follows. Given an automaton state p = ■ ■ -<S>b n ®b n+ i^)- ■ ■ , 
let c s . . . C2C1 be the subsequence of . . . b n -ib n b n+ i . . . obtained by dropping 



all the bj's such that bj 



Each Cj has the form 



= 1. 



. . lrrii . . . rrik 2 < nil < • • • < < M + 1 for some 



k > 1, 



for which we set 



mi — 1 • • • rrik — 1 



G BL 
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Let T(p) := (((ci -Ci) -03) c s ) be the semistandard tableau constructed 

from the successive products of Cj's defined via the row insertion as in |Fj] 
p. 11. By virtue of the C/'(A^_i)-invariance |[l'OY[ |, it is a conserved quantity 
under any time evolution T K , i.e., T(p) = T{T K {p)). In the context of the 
Robinson-Schensted-Knuth correspondence, T(p) stands for the P-symbol. 



For any 1 soliton state p = 
also check that T(p) equals 



(0), b e B' k , one has T{p) = b. One can 



2 2 3 



1 


1 


2 


2 


3 


3 





in Examples [2TT| , pTS| and |3.11| , respectively throughout the scattering. 



4 Automaton as Ultradiscrete KP equa- 

tion 



Here we investigate the automaton constructed in Section ^]3| ^from the 
viewpoint of ultradiscretization [TTM^ , [MSTTTfl . With the same notations 



as (|T|) we define u^j and ■ to be the multiplicities of (M + 2— j)th content 
of b l n and v n , i.e., 



y n ~~ K u n,M+l, u n,M, ' ' ' , U n,l), 



V n ~ ( V n,M+l, V n,M, ' " , U n,l> 

Proposition 4.1. T7ie map 



(21) 



expressed by (1 < j < M) 

~ U ™J = max [^l - On, X 2 -On,"' , X 3-l ~ On, Xj - K t , - • • , X M - K t , 0] 

— maxLYi — 8 n , X 2 — 9 n , ■ ■ ■ , Xj — 9 n , Xj + \ — K t , - • ■ , X M — K t , 0], (22) 

(23) 



V n+l,j U n,j + V n,j U n,j , 



M 



where Xe = X^. e : = Vii^ + ^J^rM- Actually ( \j2$) is valid also for j 
M + l. 



i=i 



i=l 



Proof. In the present proof, we abbreviate u n j and i>* ■ to and Vj respec- 
tively. We also put Uj + M+i = Uj, Vj+M+i = Vj etc., i.e. each suffix is defined 
modulo M + l. 

We define uf\ vf ] (j = 1, 2, M + 1, k = 1, 2, M + 1) as follows. 
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(1) Let Auj = Avj + i := min[u.,-, u J+ i], and := Uj—AujCvj 1 ' := Vj—Av 



,(l) - 



for j = 1,2,...,M + 1. 



(2) For Vj, we define Au^ = Auj^ 2 := min[^ iJ , u^ 2 ], and 



(i) 



.(2) 



u:- := u:- 



(i) 



A C 1 ) ( 2 ) (!) A 

Auj , Vj := uj - Auj ' . 
(3) Similarly to the step (2), we recursively define An * = Av^ + P 



,(!) 



. r (£-1) (i-lh {£) 

mm|uj ,v y j+e '\, Uj 



for t = 2, 3, 



M + l. 



u^ £ ^ — Au^ ^ and uj^ 



(£-1) a (£-1) 

uj - Auj 



From Proposition |2.3| , we see that u^ M+1 ^ and uj M+1 ' ) are the numbers of 
unconnected dots in (M + 2 — j) th box in the column diagrams for b l n and 
v\ respectively. See Example pT4] Noting that Aui = Au^- we have 



<} = «i + «} + j - «r +ij = «i + «j j - «r ( 24 ) 

for 1 < j ' < M + 1. The following formulae are easily shown by induction: 



= max 



£-i 



£-1 



Y u ^ Y u i +i + + ^ " ■ ' + v j+*> Y 



>i+i 



8=0 



-u J+ i — max 

"£-1 



i=l 

£-1 



i=2 



i=l 



8=1 



8=1 



£-1 



>j+i 



uj = max 



8=1 

£-1 



i=2 



8 = 2 



£-1 



i=2 

£-1 



(25) 



8=0 



-u 7 _i — max 



8=1 

'£-1 



8=2 



i=l 



8=1 



8=1 



£-1 



£-1 



Y v o~^ Y V i~ i + u 3~ 2 ' ' ' ' ' v 3~£+i + 



u 



J-8 



8=1 



8 = 2 



8=2 



8=2 



•J-t 



(26) 



Noticing Uj^i = Uj+M, Vj+i = Vj-M, we find 



(M) (itf) 

uj —Vj — max 



A7 



Af 



max 



Y, u i+*> Yj u i +i + v 3+i-> Yj u i +i + X v 3+ ij 

8=0 8=1 8 = 2 8=1 

Af-1 M 
^ U j+i + ^ U j+i , U j+A f + ^ U i4 

8=M— 1 8=1 8=1 

M M M 2 

_ 8 = 8=1 8 = 2 8=1 

M M-l M 

Y v i' i + Y u i~ h v i~ M + Y u i-* 

i=M-l 8 = 1 8 = 1 



(27) 
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3-1 



M+l 



Subtracting ^^Wj + t>i from the both max[...] terms in the left hand 

i=0 i=j+l 

M 

side of the equation and using the relations: um+i = d n — Uj and v m+i = 



3=1 



M 



Kt — Vj, we get (|22 ) from (|24|). Noticing that the number of dots of two 
i=i 

column diagrams are preserved in the rule, we obtain (|2"3"|). □ 
Our goal in this section is to show that (22) and (23) are ultradiscrete 
limits of the (one-constrained) nonautonomous discrete KP equation (ndKP 
eq.): 

(b n+1 - c j+1 )r(t,n,j)r(t + 1, n + 1, j + 1) + (c j+1 - a t+1 )r(t + l,n + l,j)r(t,n,j + 1) 



+(a t+ i - b n+ i)r(t,n + l,j)r(t + l,n,j + 1) = 0. 

Here arbitrary complex parameters. The ndKP eq. (f28"|), which 

is sometimes called the (nonautonomous) Hirota-Miwa equation, is equiva- 
lent to the generating formulae of the KP hierarchy H, |DJKM||. Its soliton 



(28) 



solutions, Lax operators, Darboux transformations etc. have been investi- 
gated in ||WTS|| . We set a t+ \ = 1 + 5 t and b n+ \ — 1 + j n . We also assume 
that Ci = 1, C2 = C3 = • • • = Cm+i = and 



r{t,n,j + M + l)=r{t,n,j). 



(29) 



The constraint (|29"D is an analogue of M-reduction of the KP hierarchy which 
restricts the space of transformation group of r functions to the subgroup 
generated by A$ [pjKM|| . Let 



U* ■ 



V* ■ 



r(t,n + l,j)r(t,n,j + 1) 
r(t : n,j)r(t,n + l,j + 1)' 
r(t + j + l)r(t,n,j) 
r(t + l,n,j)r(t,n,j + 1) 



(30) 



for 1 < j < M. We also introduce a small positive parameter s, and put 
5 t = exp [— K t /e] and 7„ = exp \—0 n /e\. Then we have 



Theorem 4.2. Let 



/ . 



t 



lim e log £71 
lim e log VI* ■ 
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be the ultradiscrete limits for 1 < j ' < M, and specify u l r 



n,M+l 



and v 



n,M+l 



by 



Y,f=i ut n,j = °n and Y,f=i vt n,j = K t- Then {u l nJ } and {v^} satisfy 
and ([Hj. 



Proof. We use abbreviations: Tj := r(t,n,j),Tj := r(t + 1, n,j), r n j : = 
r(t, n + 1, j), r n j := r(t + 1, n + 1, j). The ndKP eq. (|28| ) with the constraint 
(|29|) is rewritten as the following M + 1 simultaneous equations: 

(1 + 7n)ri< 2 - (1 + StjT^ + (S t - ln)r n ,i4 = 0, 
(1 + 7n)r 2 < 3 - (1 + 5 t )r* )2 T 3 + (5t - 7nK, 2 T^ = 0, 

(1 + 7n)l"MT* iM+1 - (1 + ^) r n,A/ r M+l + (S t ~ ln)r n ,M T M+l = °> 

7 n r A f + ir* )1 - 5tT* iM+1 n + (5 t - -f n )r n ,M+i T i = °- ( 31 ) 



Defining 



X l := r n,l r 2 r 3 •••Tm+1, 1/1 : = T"n,l r 2 r 3 r 4 ' ' ' Tm+X, 
%2 '■= Tl'^ i 2' r 3 •••T M+ i, 2/2 := Tl7"n,2 r 3 T 4 • • 



: = T1T2T3 •••T*j U+ i, VM+1 '■= T~It 2 T3T4- ■ -T nt M+l, 

x := (xi, x 2 , • • • , x M+ i) T , y := (yi,y 2 , ■■■ , yM+i) T , 



we obtain 
where 



Lx = (S t - 7„)y, 



/ (1 + 50 "(l + 7n) 

(1 + S t ) -(1+7, 















\ -7n 

Its inverse matrix is easily calculated as 
L 1 



(1 + St) 




D/ ((l + S t ) M S t -(l + ln ) M ln ) 








\ 



'(l+7n) 

s t J 



(i + ln ) M +^(i + s t y-\ 



j = M + l, 



hi 



S t (l + S t ) M+ ^-\l + 7n) J -\ J > i 0V M + 1), 

7n(l + 7n) M - J+J (l + StY~ J -\ j<i-l (j^M+1). 
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Thus, for < St, In <C 1, we have 





( 6t 


s t 


s t ■ 


■ s t 


1 \ 




In 


s t 


s t ■ 


■ s t 


1 


{5 t - ln )L- x ~ 


In 


In 


s t ■ 


■ s t 


1 




V In 


In 


In ■ 


■ In 





Precisely speaking, A ~ B means lim e\ogA{e) = lim e log B(e). Since 

e— >+0 



we have 



3-1 



M 



X 3 ~ ln^2yi + S t ^2 V* + VM+1, 
i=l 

3 



i=3 
M 



Xj+1 ~ ln^2yi + s t ^2 V* + VM+1, 
i=3+l 



i=l 



7n E»=l(^/l/M+l) + 5* YliLjiVi/yM+l) + 1 

rs-< ■ 

X 3 + l 7n J2i=l(Vi/ VM+l) + S t J2iLj+l(Vi/ VM+i) + 1 



(32) 
(33) 



(34) 



From the definition of U n j and V n j, we find that the left hand side of ( |3~4"D is 



M 



equal to U^/V^j and that (yj/yM+i) — Y\_ ^n,i H ^M" ^^ nce ^ holds that 

t=j i=i 



lim e log ' 



lim £ log [right hand side of (34])] 



we have (^) by putting 



lim elogC/* 



3 ' 



71, J 



From the definitions (|3~0D, we have 

rt+1 



lim e logy* -. 



y* . 

n,] n,j 



IP V 1 

n,j n+l,j 



which gives (^) in the ultradiscrete limit. 



□ 
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Next, we consider soliton solutions to the automaton. It is obvious 



that if the limit : 



lim e\ogr(t, n, j) 



exists, then from (|30| ) we have for 1 < j < M 



u 



":l 



yt _i_ yt yt yt 

n+l,j ' n j+1 n,j 



J'+l' 



v.„ 



(35) 



(36) 



. yt+1 _|_ yt _ yt+1 _ yt 

From Theorem |4.2j , they satisfy (p2|) and (^). Hence we have only to know 
3^ • to get solutions to ( ^2|) and (p3|) . We will call YjJ • an N soliton solution 

to the Ay} automaton when it is an ultradiscrete limit of one parameter (e) 
family of a certain M x N soliton solutions r(t,n,j) to the ndKP eq. (|28|) 
as explained in Appendix ||. It indeed corresponds to an N soliton state in 
the sense of Section |3~T| . 

The following fact is well known [ DJKM , WTS . 



Proposition 4.3. The N soliton solution to (j&h ) is given by the vacuum 
expectation value: 

r(t,n,j) = (vac\g(t)\vac), (37) 



N 



g{t) = Yl (1 + a k ij(p k , t)i[>*(q k , t)) . 



(3? 



k=i 



Here t = (t, n,j) and a k (k = 1, 2, • • • , N) are arbitrary complex constants. 



ip{p,t) 



i'=i 

3 



. t' n 1 j'=l 



with 



]ln'=l -^n' 

1 



1 < n 
n = 

nL +1 ^ ! n<-l, 

and if) (p), ip*(q) are fermionic field operators which satisfy 
{ip(p),ip(p')} + := ip(p)ip(p') +ip(p')ip(p) = 0, 
{^(g),^(g')}+ = 0, {V(p)^*(?)}+ = for(p^g), 

(t»ac|'0(pi)'0(P2) ■ ■ ■'4 ) {pr)'4>*{qr)'4>*{qr-i) ■ ■ ■ \vac) 



det 



Pi Qj / l<i,j<r 



Ui<j(Pi - Pj)(Qj - g 
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The iV soliton solution (|38| ) is also a solution to ( |3~1~D when it satisfies the 
constraint (|29|). We can easily show 



Proposition 4.4. T7ie constraint is achieved if it holds that 

= 1 (A; = 1,2, •• - ,N). (39) 



\ M / 1 



PkJ \l-Pk 

Note that, for a given p*., there are M q^s which satisfy fl3"9] ) and q% 7^ p^. 
We use this fact to construct explicit solutions. 

From Propositions I4.3I and POL we can construct a class of iV soliton 

Theorem 4.5. 



solutions to the A M automaton. The result is summarized as 



N 



(40) 



J2^K (l) (t,n l3 )~A(fI;j) 
_«=i 

is an N soliton solution to the automaton. Here fl = (/ii, fi2, Hn) 
(fii = 0,1) and max[- • • ] denotes the maximum among the 2 N values obtained 

by putting fii = or 1 for % = 1, 2, • • • , N. 

j t n 

K®(t,n,j) = K$> ~ X>n[«*,L«] +E mi *'> L(i) ]> 

j'=l V n' 

where the sums here are generally defined by 

n [ YuZ=l X n' 1 <n 

^X n , := I n = 

I -J2°n'= n +l X n' U < -1. 

ij 1 ^ (1 < i < N, 1 < j < M) are non-negative integers which satisfy 

M 

i=i 

L« > > • • ■ > 

£f>ef>...>£f\ ( J = 1,2,---,M), 
and Kq^ is an arbitrary integer. In the case: ' ^ l ' * *i>*2j 'V 



i/ie phase factor A(jl;j) is given by 



fii = otherwise, 



j) = J> - 1)L<**> + J2 ( X(tk) (j +k-l)- X<M(j)) 



k=l k=l 
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wh 



7« 



with 



/>(*) - />(*) 

l j+M - l j ■ 



The proof of this theorem is parallel to that in [ TTM]. We give the detail 
in Appendix [FJ. For AT = 1 it is the general solution, and we conjecture that 
it is also so for M = 1. Except these cases the above result does not cover 
the arbitrary initial condition. There is some freedom to employ different 
'phase factor' A(p,\ j) than the above one depending on the way in taking the 
ultradiscrete limit. 



5 Summary 

In this paper we have introduced the A$ automaton, which is a crystal 
theoretic formulation of the generalized box-ball systems. In terms of the 
box-ball systems, it corresponds to the dynamics of M kinds of balls, where 
the carriers and boxes have arbitrary and inhomogeneous capacities. We have 
introduced the solitons labeled with the crystals B' k of U'{A^_^). Scatter- 
ing matrices of two solitons are identified with the combinatorial R matrices 
of Z7'(A^_ 1 )-crystals. Piecewise linear evolution equations are obtained and 
identified with an ultradiscrete limit of the nonautonomous discrete KP equa- 
tion. It allowed us to construct a class of iV soliton solutions. We have left 
the studies of phase shifts in the scattering and construction of iV soliton 
solutions corresponding to arbitrary initial conditions for iV > 2 as future 
problems. The interplay between the ultradiscrete limit of the classical in- 
tegrable systems and the q — > limit of the quantum integrable systems 
elucidated in this paper deserves further investigation. 

Acknowledgements. The authors thank M. Okado and Y. Yamada for 
discussions and sending a preprint prior to the publication. Thanks are also 
due to J. Matsukidaira, A. Nagai, J. Satsuma and D. Takahashi for helpful 
discussions about soliton solutions. 



A Proof of Proposition |StB 



First we show that it suffices to prove Proposition |3.8| for M = 1 and h = k. 
Without a loss of generality we may set M = 2 and consider the time evo- 
lution T K=OQ . We find it convenient to adopt the equivalent box-ball system 
picture explained in Section [2.4| . Thus the elements in Bf c in ( |13"D will be 
represented as . . . 131. .2 . . . for example. It stands for the array of the balls 
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with the indices 1,3,1 and 2 and . denotes an empty box. (So they do not 
correspond to the letters in the semistandard tableaux in the crystal nota- 
tion.) We keep the same notation 9 to denote the map corresponding to fll^D 
in the box-ball picture. It groups the array of balls and empty boxes locally 
together into the boxes with capacities ... , 9 n , 8 n +i, .... Then the assertion 



of Proposition 3.8 is that the scattering 



k—h l+h—k 




(41) 

takes place for sufficiently large t. Here Ti,T% are the ball- moving operators 
defined in Section [2J], and we have used = (T 2 Ti)* in view of (|i~0f ) and 



the fact that the balls with index > 3 are absent. In ([II]) the sequences .... 
of the empty boxes are sufficiently long since both sides are to represent the 



asymptotic 2 soliton states in the sense of Section |3.1| . Now we make use of 
the relation ^Ti) 4 = T2{T\T2) t ~ 1 Ti. From the definition of the operators 
Tj's and the assumption that the 2 solitons are enough separated, (f4"T|) is 
equivalent to 

k l+h—k k—h 

§(. . .2~^~2. . . .2---21---1. . .). 

(42) 




But this is justified once one establishes 

l k 




(43) 



because ([12]) and (|43| ) correspond to the same canonical system 



9(. . . 12 • • • I I + 1 • • • I + k . . . ) in the sense of Section |2]4] with respect to 

the relevant time evolutions and therefore they possess the parallel time 
evolution pattern owing to (|i~l~l). In this way the proof of Proposition |3.8| is 
reduced to (|43|), which is equivalent to the case M = 1 and h = k. 
Now setting LP^ = I and = k, we are to show 

Proposition A.l. Set M = 1, assume that n t ^> (it) and 9 n < 
for all but finitely many n's. Then two solitons with amplitudes and 

L m (L (d > L ( 2 ) j 

scatter into two solitons with amplitudes and 

respectively. 

Namely, the amplitudes of two solitons do not change after the collision. 
To prove the proposition, we need several lemmas. The following two lemmas 
are obvious. 
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Lemma A. 2. For given integers K\ and K 2 , if there exists an integer n 
such that 



no n 



Ki + Y^ min [6 n ,, > > K 2 + min I *, L {2) ] , 

n' n' 

then, for n > n , 

n n 

K. + Y, min [6 n ,, >K 2 + ^ min [6 n ,, , 

n' n' 

and for n < no 

n 

> if 2 + ^min [6 n ,,LW] . 

n' 

Lemma A. 3. For given integers K\ and K' 2 , if there exists an integer n 
such that 

no n 

K' 2 + Y mm [9 n >, L (2) ] > 2L (2) >0>K\ + J2 min [9 n ., L (1) ] , 

n' n' 

then, for n > n , 

n n n 

K\+K' 2 +J2^ [^,L (1) ]+^min [0 n >,LV>]-2LV> > i^+^Tmin [6 n ,,L^] , 

n' n' n' 

and for n < n , 

n 

> K\ + ^min [0 n /,L (1) ] . 

n> 

Now we define an integer N (t) for given integers K 2 and t as 

N (t) 
n> 

> 

N (t)-1 

>K 2 -L ( - 2 k+ J2 min[0„,,L( 2 )]. 

n> 

With this N (t) we can show 
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Lemma A. 4. For any integers K\ and K 2 , we have 

N (T) 

lim (K 2 - L^T + V min \B n >, L (2) l 



n' 

-Ki + L (1) T - min [9 n ,, L (1) ] ) = +00. 

n> 

Proof. From the definition of Nq(€), we have 

N (t) 

-L&> < -tL® + Yl min P»'> L(2) ] < L(2) - 

n'=7V (0)+l 

Hence we have 

N (t) 

A(t):=-L ( - 2 k + Yl min[^,L( 2 )] 

n'=AT (0)+l 
/ Wb(t) 

_ _ L (D t+ ^ min[^,L«] 

y n'=7V (0)+l 
N (t) 

= t (L« - L^) - ( min - min l>« 

n'=Af (0)+l 

> t (L« - L®) - £ [O n ,-rmn[e n/ ,LW]) 

n'=7V (0)+l 
N (t) 

> tL (D _ _ g ^ 

n'=7V (0)+l 

From (0), we obtain an inequality: 

* >_1 + 7W min[^,L (2) ] . 



Thus, from (1451), we find 



LP) 

n'=JV (0)+l 



A(t) > -(L (1) + L (2) ) + ^ min 

re'=JV (0)+l 
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LW - LP) 
L( 2 ) 



^,LW 



Since L (1) > L (2) , L (1) > 6 n for all but finitely many n's and lim No(t) = 
+00 which is seen from (fj4|), we find 

lim A(t) = +00. (47) 

This suffices to prove the lemma. □ 
Now we prove Proposition [A.l| . From (^) we have 

?/ ■- v t+1 

u n • u n,j=l 



Vt+1 _ yt+1 _ v t+l 1 yt+1 ( A c 
I n+l,l I n+1,2 I n,l ' 1 n,2 • 



Specializing Theorem [15] to a two soliton solution with M — 1 and Kt = +00, 
we have 

Yn+1,1 = max [0, K x (n, t), K 2 (n, t), K^n, t) + K 2 (n, t) - , 

K n *+J )2 = max [0, Ki(n, t) - L {1 \K 2 (n, t) - L (2 \ K^n, t) + K 2 (n, t) - - 3L (2) ] 

n 

Ki(n, t) = Ki + J2 min Pn', - ^ W (i = 1, 2). 



n' 

"t+1 



Note that Y^ 2 = due to the last equation in (^T|) and the condition 

8 t = exp[— K t /e] = 0. Given L^ 1 ' > L^, there exist integers n x , n 2 , j, r 1; r 2 
that satisfy 

ni<ni+j<n 2 , 1 < n < min(L (1) ,# ni ), 1 < r 2 < min(L (2) , 6 n2 ), 

ni n\ 
Kt + ]T mil1 [ n>,L {1) ] > > K 2 + Y min Pn', ^ (2) ] , 



- + ]T min [0 ft ,, L«] > > K 2 - + ^ min [0 n ,, L< 2 )] , 

n' n' 

where 1^ (i = 1,2) is defined by 

K x = t x - J] min p„/,L (1) ], 

n' 

ir 2 = r 2 + 2L( 2 )-^min [^,L«]. 

n' 

From Lemma |A.2j , we find at £ = that 

F^i,! = max [0, K x {n, 0), ^(n, 0) + K 2 (n, 0) - 2L^] , 

y^ lt2 = max [0, Ktin, 0) - L«, K x {n, 0) + fC 2 (n, 0) - - 3L< 2 >] .(49) 
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Substituting ( pf9| ) into eq. P8|), we obtain 



On 



Tp) - V n ' 1_1 



n 



r 2 

0„ 



(2) _ yn' r l g _ 

Z-~in —112+1 n z 



for n < n\ 

for n = ni 

for ri\ < n < n\ 

for n = n'\ 

for n'i < n < n 2 

for n = n 2 

for n 2 < n < n' 2 

for n = n' 2 

for n' 2 < n, 



(50) 



where n'i (i = 1, 2) are defined by n\ = rii + 1 if = Lw, and otherwise by 

n'i ri'j— l 

L {1) - Yl n >-ri<O<L®- d n'~n. 

n'=7ii+l n'=rii+l 

Thus we see that the two soliton solution can correspond to any initial con- 
figuration in which L^- 1 ' soliton is situated left hand side of soliton with 
sufficient spacing. Hence, to prove the proposition, we have only to show 
that the solution u l n describes the two soliton state in which soliton is 
left hand side of l/ 1 ) soliton for t ^> 1. 

From the definition of N (t) and Lemma |A.4| , there exists T and j such 
that 

N (T) JVo(T) 

K 2 - L {2) T + Y min [0*', L{2) ] >0>K X - L (1) T + ^ min L {1) ] , 

n' n' 

N (T)+j N (T)+j 

K 2 _ L m T _ L m + J2 mm[e nl ,L^]>0>K 1 -L^T-L^+ ]T mi n [0 n ,,L«] 



Thus, from Lemma A.3 , we have at t = T that 

= max [0, K 2 (n, T), i^(n, T) + K 2 (n, T) - , 

= max [0, K 2 (n, T) - , K^n, T) + if 2 (n, T) - L« - 3L< 2 >] . 

Substituting these into eq. (pESD, we find that describes a configuration in 
which soliton locates around n = Nq(T) and soliton does around 
n ^> N (T). This completes the proof. 
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B Derivation of N soliton solutions 



Here we explain the derivation of the N soliton solution in Theorem |4.5| along 
the simple cases N = 1 and N = 2. First we consider one soliton solution. 
We will show that it has the form: 



y*. 



max 



i-i 



t-i 



n-l 



mm[K t ',L\ + ^2 min[^ n /, L\ 



i=l 



(51) 



where L is the amplitude, K is an integer which is related to the phase 
of the soliton, and (i = 1, 2, • • • , M ) are the non-negative integers which 



M 



correspond to the number of zth balls in the soliton and Vj^j = L. We 



»=i 



give some details of its derivation, because similar technical difficulties in 
obtaining multi-soliton solutions are resolved in the same way. 
To obtain (pT]), we take g{t) in fl3"gp as 



M-l 



g(t) = H (l + c e (p)^(p,t)r(qi,t)) 



E=0 



l + lP{ P ,t)(f>*{p,t) 



(52) 



M-l 



1=0 



where qe (£ — 0, 1, - • • , M — 1) are the roots of the algebraic equation : 
x M (l -x) -p M {l-p) 



x — p 



0, (x ^ p) 



(53) 



for a given real number p [(1 + M -1 ) -1 < p < 1], and q(p) (0 < £ < M — 1) 
are complex coefficients which will be determined later. Since ( |53"D has one 
real positive root, we assume that go is positive and we put r] = qo/p. Then 
p and go satisfy 



V 



M 



1 — p = T] 

% = 1] 



M 



1 — T] 



1—7] 



rj 1 



7] 



M+l 



(54) 
(55) 
(56) 
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The r-function r(t,n,j) is given by vacuum expectation value as 
r(t,n,j) = (vac\g(t)\vac) 

_ ^ „ i ^y-^ 1 / i- P /(i + * n A-W(i + 7o) \ 

h P-U\PJ U-^/(l + S )J U-p/(l + 7o)J 

H / l-p/(l+^) \ ^ A-g//(l+7nQ \ 

11 l- a,/(i + s#) M 1 l- p/fi + Tn'i ; ■ 1 J 



X 

f 



We introduce a small positive parameter e and put rj = exp[— L/{Me)]. We 
also put 



ceip) fqtY™- 1 fl- P /(l + 5 ) \ /l-%/(l + 7o) 



Q(P): P~qe\pJ \l-q e /(l + 5 )J Vl-p/(l + 7o) 

JVo 

x n n (i-w(i+7n0), (ss) 

t'=-T n'=l 
M-l . x s 

E ^) ( f ) ( SGZ )' ( 59 ) 

£=0 ^ P ' 

where T = T (e) and N = N (e) are positive integers which satisfy T ~ 
iVn — 1/e. Hence, lim Tn = lim Nn = +oo. Since 

M-l / N s+M 



M-l 



M-l / x s / , 



v V i - % 



£=0 

oo M— 1 y- v s+i 

i=0 £=0 ^ ' 

oo 



we have 



M— 1 / oo 



Xp (s + M) = X> - P) m P*"ft(0 + 0, (60) 



i=0 \£=0 
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where Qo(i) — 1, Qi(i) — i + 1 and 

(e-i)M+i ki fc^-i 

«w= E E "E 1 

k!=(e-l)M k 2 =(£-2)M k e =0 

=^n(^+<+j+i), 

for £ > 2. Note that Xp( s ) is a re& l function when Xp(j) (0 < j < M — 1) are 
real. The ratio Qe+i(i) / Qe(i) {£ > 1, < i < M — 1) is calculated as 

(£ + l)(M + l)+i-A- / M 



fc=2 

<(M + 1) (l + I 
< (M+ l)e. 

Hence, if it holds that (1 - p)p M < (M + l) -1 e _1 , we obtain 

M-l 



| Xp(s + M )| <(1-P)E( 1 + « + ^ i-ft-p^ + De ) !*<« + 

(61) 



M-l 



Thus we find x p (s + M) ~ Xp( s + f° r sufficiently small rj. 

We assume the following for XpU)'- 

XpiX) = Xo 

Xp (2) = N iy ^ Xp (l) 

Xp (3) = N^XpW 

Xp (M) = N M _ iy ^x P (M - 1). (62) 

Here xo is a positive number which is related to the initial phase of soliton, 
y = exp[— 1/e], £j and Nj = Nj(e) (j = 1, 2, • • • , M — 1) are non-negative 
integers and positive numbers respectively. They are also supposed to satisfy 

M-l 

£ M := L - t-i > 0, 

3=1 

lime log iVj-(e) = 0, (63) 
£ N^<e N \ 



38 



for a sufficiently large positive integer N* . From these conditions, ce(p) 
(0 < I < M — 1 ) are uniquely determined by the equation: 



/ 



9o qi 

9o 9? 



\9 M 



9m -i 



9M-! / 



/ 5>(p) \ 
\ c M -lip) ) 



( PXpiV \ 
P 2 X P (2) 

^ p M x P iM) ) 



(64) 



Note that the determinant of the M x M matrix in the left hand side is equal 



to 



From 



■M-l \ / \ 

IT ) ( n^-? — 9*) ) ^ ^" ^ snomc i be also noted from (^xL 

. i=0 / \j>i J 

Xpif)>e- N \ P ii + l) forVz, 
Xp(i) > Cexp[L/e}x P ii + M) for Vz and 3C > 0. 

we have 

M-l 



that 



(65) 



n-:i 

X 



j t-1 



f=-Tn 



-1 t-1 



P 



+ 4 



n(i-rfe) II - T 



-1 Wo 



+ 7r 



-1 



(66) 



Hereafter we restrict ourselves to the region: |n| < iVo and |t| < T . Noticing 
that 



t-i 



t'=-T 
= 1 



n i-ifd n I- 



-1 JVo 



n'=n 



+ Jr. 



( 1 + 5, J + S ( 1 + 



=: 1 + ai 



t=-T 
9f 



P 



<5* 
+ a 2 



n—n 
2 



+ a 3 



+ 7n< 

'ft 
,P 



ft + • 



we find 



t— 1 / \ n— 1 /■ 



p 



+ 7r, 



+ 



-1 oo 



J^OiXpO' + Oj ( 67 ) 



i=0 



oo 



where ao = 1 and a i+ i/ai ~ e \ From (|65|), we have < ^^aiXp(j + *) < 

i=l 

XpU) f° r sufficiently small e. Putting xo = expfffo/s] and noticing the 
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relation: 



lim elogfl — p) — —L, 

3-1 

Ifin elogx P 0') = K - 

i=l 

( v V 1 

lim e log 1 = min \L, 9 n ] 

lim £ log ( 1 - — — ) = - min [L, K t \ 



we obtain 



lim e \ogr(t, n, j) = max 

£— >+0 



j-l t-1 n-1 

min^j/, L] + ^ min^', L] 



(6* 



Since lim AT (e) = lim T (e) = +oo, we have shown that (pll) is a one 

soliton solution to the A$ automaton. 

Next we consider two soliton solutions. From the above arguments about 
one soliton solution, we see that the field operators ip(p) and (f)*(p) are essen- 
tially determined by L, ij (j = 1,2, •■ ■ , M) and K . Therefore we denote 
these operators by 

^(jp)=^{L:e), <P(jp) = <i>*(L;{e j };K :e). (69) 

Then we take 

g(t) = (l + V(pi,*)0*(pi,t))(l + V'(P2,t)0*(P2,*)), (70) 

where 

=1>{L®:e), <F{p i ) = <P{L^i®y,K®:e) (i = l,2). (71) 

We also assume > L (2) and if > if (j = 1, 2, • ■ ■ , M). As we shall see 
below, the latter condition turns out to be a natural constraint for soliton 
solutions. Using the similar notations as above, we have 

r(t,n,j) = («ac|(l + ^(pi,t)0*(pi,t))(l + ^(p2,*)0*(p2,*))|vac) 

= 1 + (vac\i/j(p 1 ,t)(f)*(p 1 ,t)\vac} + (vac\ip(p2, t)0*(p 2 , t)\vac) 
+(vac\i/>(p 1 , t)<f>*(p!, t)V»(p2, *)0*(P2, t)|uoc). (72) 
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The second and third terms are calculated in the same way as above. The 
fourth term is evaluated as 



(vac\ip(pi, t)(f)*(pi,t)il;(p2, t)<p*(p 2 , t)\vac) 

M - lM - 1 ( \ pi - p 2 )(q (2) •■ |J - 1 - 



/Ji)\ j t-1 / (i) \ ~ l t-l , 

n-l / \ -1 JVo / (t) ^ _1 



X 



n' n=i 

We define by 



M-l / (i)\ s 

X«W-E^) 7 (< = 1,2), (74) 



£=0 



and suppose 



X P M = N?y^ XlH {2) 



Xpi (M) = nS^m-ixpAM - 1), (75) 

where positive numbers satisfy the similar inequalities to flB3|). From 
the assumption: $p > £^ (j = 1, 2, • • • , M), it is always possible to choose 
Nf ] such that 

XpJJ + jj > x P i(j + 1) / 76 x 

Then ([73]) is expanded as 

2 t-1 n-l / \ / \ -1 



a-^nnn('-i^)('-i^)' 

oo oo 

X] ( a ^'Xpi(j + i)Xpa(j + 1 + - Ki'XpS + i)Xpi(j + 1 + 0) > 



X 

i=0 i'=0 
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where the coefficients are defined by 



V^V\ \ TT TT / 1 q lk \ TT / i q l 



t-i / Jk) \ - 1 No / (k) 



and fc^j' = a «',i- F rom (|65f) i we evaluate 



oo oo 



a 0,0Xpi 



i=0 i' =0 

oo oo 



& o,oXp 2 (j)Xpi(j + 1) > XI hi'X P2 U + i)Xtn(J + 1 + 0- 

i=0 j'=o 
i+i'^0 

Then, noticing a 0) o = P2> ^o,o — Pi an d using (|76|), we find 

lim elogr(* + l,n + l, j + 1) = max [0, K {1) (t, n, j), K {2) (t, n, j), 

K^(t,n,j)+K^(t,n,j)-A(j)], (77) 

K®(t,n,j) := K® -j^if - ^min 

i'=i *' 

n 

+ 53min[0 n /,L«] (i = l,2), (78) 

re' 

A(j):=L (2) +4?i (0<J<M-1). (79) 

This gives a two soliton solution. For the scattering where the larger soliton 
overtakes the smaller one like (I) in Theorem \5.10[ the integer £j (1 < j < 
M) corresponds to the number of j th balls in the larger soliton at t — > — oo, 

(2) 

and £j corresponds to that of the smaller soliton at t — > +oo. Since the balls 
in the smaller soliton at t — > +oo must be included in the larger soliton at 
t —y — oo, the condition tp > €p must hold for soliton solutions. Similarly, 
for the scattering where the smaller soliton overtakes the larger one like (II) 
in Theorem |3.10| , the integer £p (1 < j < M) corresponds to the number 

of j th balls in the smaller soliton at t — > — oo, and £j corresponds to 
that of the larger soliton at t — * +oo. We should also note that there are 
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several freedoms to choose the 'phase' A(j) in taking the ultradiscrete limit. 
However we conjecture that the above choice will cover all the canonical 
systems, hence essentially all the time development patterns for N = 2. 

The N soliton solution flUf ) is obtained in the same way. The key in the 
construction is to evaluate the expansion: 

(vac\ip(p 1 )ip*(q 1 )ip(p2)ip*{q2) • • • ip{pr)ip*{q r )\vac) 

= (vCLC^ip^ipz) ■ ■ -^{PrW{qrW{qr-l) ' ' ' ^*(<?r) \vac) 
_ ril<i<i<rbi-^i)fe -^) 



Ul<i,j< r (Pi ~ Qj) 

(<?r _1, ?r-i ' ' ' 92 + other terms) 



UlKiKjKriPi ~ Pj) 



YYi=lPi 

and show that this term gives the phase factor A(jl;j) and the "other terms" 
do not contribute to the final results. This can be done in the same manner 
as in the case of two soliton solutions. We take 

N 

9(t) = YlO-+i>(Pi,t)<f>*M), (80) 

i=i 

where 

^(R)=^(L«:e), ^( Pl )=^(L^;{£f};K^ : e) (i = 1, 2, • • • , N). 



ID 



We suppose 
and 



L« > > . . . > iW, 



Zf>lf>--->lf\ (j = l,2,---,M). 

Note that this implies: p\ > P2 > ■ ■ ■ > Pn- The latter condition is also 
a natural constraint for iV soliton solutions as in the case of two soliton 
solutions. Finally we find that the result is given by (f40l) . 
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